Mapping the X-Shaped Milky Way Bulge 
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ABSTRACT 

We analyzed the distribution of the RC stars throughout Galactic bulge using 2MASS data. 
We mapped the position of the red clump in 1 sq. deg. size fields within the area |^| < 8.5° 
and 3.5° < \b\ < 8.5°, for a total of 170 sq. deg. The red clump seen single in the central area 
splits into two components at high Galactic longitudes in both hemispheres, produced by two 
structures at different distances along the same line of sight. The X-shape is clearly visible in 
the Z-X plane for longitudes close to ^ = 0° axis. Grude measurements of the space densities of 
RC stars in the bright and faint RC populations are consistent with the adopted RC distances, 
providing further supporting evidence that the X-structure is real, and that there is approximate 
front-back symmetry in our bulge fields. We conclude that the Milky Way bulge has an X-shaped 
structure within |^| < 2°, seen almost edge on with respect to the line of sight. Additional deep 
NIR photometry extending into the innermost bulge regions combined with spectroscopic data 
is needed in order to discriminate among the different possibilities that can cause the observed 
X-shaped structure. 

Subject headings: Galaxy: bulge — Galaxy: structure — stars: distances — stars: late-type 



1. Introduction 

In the last decade consensus has been reached 
that the Galactic bulge has a boxy s hape in the 



COB E/DIRBE near infrared map (JDwek et al. 



However, recent analysis of the RC stars across 
the bulge a rea suggested a more complex X-shaped 
structure dMcWilliam et all [2OI0I: Izoccalil I2OIOI: 



I995I ) because it contains a bar. This bar is tilted 
by ^ 20 — 25 degrees from the Sun-Galactic center 
direction, with axis ratio s 1 : 0.33 : 0.23, and a, scale 
length of ~ 1.5 kpc (e g. IStanek et al."l994' Il997t 



.JStan 
12005; 



Babusiaux fc Gilmord 12005; Rattenburv et al. 
20071 and references therein) . Most of these stud- 
ies make use of horizontal branch red clump (RC) 
stars as tracers of the bulge 3D structure, due to 
their narrow range in color, brightness, and well 
calibrated absolute magnitude. 



Nataf et all l2010t iMcWilliam fc Zoccalil \20ldi ). 
Along to / = 0° axis, both above and below 
the plane, the RC appears to split in two com- 
ponents. Both components become closer and 
finally merge at the latitude of Baade's Window 
(6 = —4°) and below, while they open up and di- 
verge at higher latitud es both above and b elow the 
plane. The analysis by lNataf et al.l poill ) is based 
on the RC dis tribution from the OGLE-II pho- 
tometry, while iMcWilliam fc Zoccalil (|2010t ) ana- 
lyzed the color magnitude diagrams (CMDs) from 
2MASS on a roughly symmetric field grid across 



the btil ge. For the field at (I b) = (0° -6°), in par- 
ticular, iMcWiniam^Zoccaii (|201o( ) had deepe 



per 

SOFI@NTT infrared photometry (Ks ~ 18 mag, 
in comparison to Kg ~ 14.3 mag from 2MASS) as 
well as deep optical data from WFI@ESO/MPG 
telescope {V ~ 23.0 mag, / ~ 22.5 mag). 

Here we complete the analysis of the RC dis- 
tribution across a much wider bulge area using 
2MASS data. The RC was mapped at positive and 
negative latitudes in 170 fields of 1 sq. deg. size, 
within \l\ < 8.5° and 3.5° < |6| < 8.5°. Our results 
are in agreement with previous studies, and ex- 
tend their findings, showing the presence of an X- 
shaped structure across the whole Galactic bulge. 

2. Observational Data 
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Fig. 1. — The bulge area analyzed here, divided in 



170 sq uare fields, overplotted on the Sc hlegel et al 



(|l998l ) dust map. The inner area was not used here 
because the 2MASS photometry does not reach 
the RC magnitude. 

This analysis makes use of data from the Two 
Micr on All Sky Survey (2M ASS) point source cat- 
alog (jSkrutskie et al.ll2006l ) for the Galactic bulge 
region within coordinates |^| < 8.5° and 3.5° < 
|5| < 8.5°. In order to map the RC position across 
the bulge, the whole area was divided in square 
fields of 1° size, centered on each integer Galac- 
tic coordinate, for a total of 17 x 10 = 170 fields 
(see Fig.[T|). Unfortunately, the incompleteness of 
2MASS in the region close to the Galactic plane 
(|6| < 3.5°) becomes rather severe at the RC mag- 
nitude, so this region has been excluded from the 
present analysis. 
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Fig. 2. — Top panel: Zoom up of the RCs in 
the 2MASS CMD of a selected field at {I, b) = 
(0°, —7°). A vertical dashed line marks the color 
cut applied in order to exclude from the analy- 
sis the foreground disk main sequence stars on 
the blue side. Contours help the eye identify the 
two clumps, also marked by horizontal thick lines. 
Bottom panel: Kg luminosity function of the bulge 
RGB near the RC, both in 2MASS and in the VVV 
data (see text). 



As discussed by iMcWilliam k Zoccalil (|2010f ). 
the most plausible explanation for the double RC 
is the presence of two structures along the line of 
sight. We therefore mapped the position of the 
RC(s) across a large area, in order to deproject 
the bulge 3D structure. The procedure was the 
following: 

i) A Kg vs. (J — Ks) CMD was constructed 
for each of the 170 fields (Fig. [2|) and a color cut 
applied, fine tuned for each CMD in order to sep- 
arate the bulge red giants (RGB) from the main 
sequence of disk foreground stars. The luminosity 
function of bulge RGB stars was constructed in 
order to identify the RC(s). 

ii) Using the RC(s) as distance indicator, we 



assume that the two peaks seen at the bottom of 
Fig. [2] are due to the presence of two overdensi- 
ties along the hne of sight. In order to map the 
overdensities in 3D space, we associate a distance 
to each of the two peaks, and further assume that 
the width of each peak is entirely due to distance 
spread. This latter assumption lead us to overesti- 
mate the distance spread, because other effects can 
contribute to the RC width. The photometric er- 
rors at the magnitude of the RCs {Kg < 13.5 mag) 
produce uncertainties of Ad ^0.1 kpc in the RC 
width, while the differential reddening should ac- 
count for Ad ~ 0.2 kpc. In addition, the RC has 
an intrinsic width due to the broad metallicity dis- 
tribution, and a possible spread in age and helium 
content, both hard to quantify. The assumption 
of the width of the RC peaks being due only to 
distance allows us to directly convert the luminos- 
ity function into a density map along the line of 
sight, by assigning a distance to the Kg magnitude 
of each star close to the RC. In this way, the po- 
sition of the main structures in the map will be 
correcu, but their width along the line of sight 
will be overestimated. The distance associated to 
each RC star was calculated as: 



M 



-5 -I- 5log{d) ^K ~Mk-A^ 



(1) 



where the absolute magnitude Afx was taken 
from the horizontal branch in the isochrones of 



Pietrinferni et al.l (|2004l ) fo r solar metallic ity, and 
Ak taken from the maps of ISchlegel et al.l ([l998). 
This makes the reasonable assumption that all 
the extinction happens between us and the bulge, 
with no dust being present between the two bulge 
RCs. 

Hi) Distance histograms, corresponding to den- 
sity maps along the line of sight, were constructed 
for RC stars in each of the 170 fields, using an aver- 
aged shifted histogram, for a bin width of 0.2 kpc. 

Figure [2] also shows the luminosity function of 
RC stars, in the same bulge field, from the Vista 
Variables in the Via Lactea ESO Public Survey 
(VVV; Minniti et al. 2010). For this figure, the 



^The intrinsic width of the RC, combined with the fact that 
we did not subtract the underlying RGB luminosity func- 
tion, could actually displace slightly the estimated position 
of the RC peak. However, this would be negligible com- 
pared to the effect of a small difference in age or helium 
content, and in any case none of these effects would change 
our conclusion about the global bulge shape. 



CASU vl.l version of the catalog has been useqj. 
In this field, the VW survey reaches down to 
Ks ^ 17.5 mag and the photometric error at the 
level of the RC is much smaller than in 2MASS. 
The presence of a double RC is fully confirmed by 
VVV data. The faint RC is stronger in the VVV 
data, due to incompleteness of 2MASS. In fact, 
the discrepancy between the two surveys increases 
at fainter magnitudes, due to confusion becoming 
stronger and stronger in 2MASS. The VVV survey 
data will allow us to map the bulge 3D structure 
down to the Galactic plane. Extensive analysis 
of the upper RGB and RC morphology across the 
bulge area covered by VVV is ongoing and will be 
the subject of a dedicated paper (Gonzalez et al. 
2011, in preparation). 

3. The X-Shaped Structure 

In what follows we will discuss different sections 
of the bulge density maps constructed above. 

The panels in Figure [3] show sections of the den- 
sity map roughly parallel to the Galactic plane, 
seen from above, each at a given Galactic latitude. 
The region within \b\ < 3.5° is excluded from the 
analysis, so the innermost sections are at b = +A° 
and b — —4°, in the middle row panels. A cross 
marks the Galactic Center, assuming a distance 
of 8 kpc (e.g., Ro = 8.00 ± 0.6 kpc; IChez et al 



20081) . The Sun would be at {l,X) = (0°,0), be- 
low and outside the limits of the plots. The map 
in each panel is constructed merging 17 vertical 
color strips for the fully sampled sub-panels, each 
showing the RC star density along a given line of 
sight, at a given longitude (x-axis). Contour lines 
mark equal density levels, in steps of 5% from the 
peak intensity. The section at & = -1-4° has a ver- 
tical, black strip due to a dust lane that prevents 
us to reach the RC in 2MASS, in this part. Due 
to the opening of the cone of view, we see nearby 
things more expanded in the sky, than far away 
ones. In order to show the distortion introduced 
by this effect, we also show diagonal thin white 
lines at constant Y linear coordinate. 

The color scale has been normalized in each 
panel so that the peak density of that section is 
yellow. This stretching of the scale allows us to 



^This is the catalog, produced by the Cambridge Astronom- 
ical Survey Unit (http://casu.ast.cam.ac.uk/), correspond- 
ing to the 1"* VVV data release, on June 2011. 




Fig. 3. — Density maps showing the structures traced by the RC near the Galactic plane, i.e., as seen 
from above, in slices of different latitudes (see labels). Individual lines of sight, at a given longitude, are 
represented by vertical strips, which are then merged together to form each panel. The central strip in each 
panel corresponds to Z = 0° with a cross marking the Galactic Center (assuming Rq = 8 kpc). The panel 
at 6 = —4° also shows the Galactic bar as traced by Rattenbury et al. (2007; white dots) fitting the RC in 
OGLE II data. The label at the bottom of each panel lists the peak value of the density histogram in that 
particular section of the 3D map. Contour plots may help the eye in regions of low density contrast. Thin 
white lines are lines of constant Y coordinate. 



trace the shape of the bulge even far away from 
the plane, where stellar densities are much lower, 
but might give the wrong impression that the to- 
tal number of bulge stars at |5| = 4° is the same 
as it is at 1 6| = 8°, which is clearly not the case. 
To clarify this, the peak value of the density his- 
togram (corresponding to the lightest yellow) is 
listed at the bottom of each panel. 

Two overdensities can be seen in each panel: 
one closer to the Sun, at positive Galactic longi- 
tudes, and a more distant one, at negative longi- 
tudes. Moving away from / = 0° to larger lon- 
gitudes, one of the two structures gradually dis- 
appear, leaving just the bright one at positive 



longitudes, and the faint one at negative longi- 
tudes. The two overdensities get closer to each 
other for sections closer to the Galactic plane, and 
they almost completely merge already at 5 = —4° 
(Baade's Window). This is the section that has 
been studied extensively in the past, and the struc- 
ture that one sees here is consistent with a single, 
elongated overdensity, i.e., the bar. Clearly, only 
for lines of sight along to ^ = 0° axis, i.e, in the 
central, vertical strip of each section, and only for 
|6| > 4° we intercept both overdensities, resulting 
in a double RC. 

Interestingly, in this map the far side of the X 
fades faster than the near side, when moving away 
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Fig. 4. — Density maps showing the structures traced by the RC in the {X, b) plane. Each panel corresponds 
to a different longitude (see labels), with ^ = 0° in the central, middle one. A cross marks the Galactic 
Center (assuming Rq = 8 kpc), with the Sun far in the left side, outside the figure, at {X,Z) = (0,0). 
Individual lines of sight, at different latitudes, correspond to horizontal color strips, merged together to form 
the panels. The color scale has been normalized so that the histogram peak in each horizontal strip has 
density=l. Note that the region close to the Galactic plane, for which we have no data, has been compressed 
here, and shown as a single, horizontal black strip. Thin white lines are lines of constant Z coordinate. The 
X-shape is clearly visible for longitudes \l\ < 1° (middle row panels). 



from the Galactic plane. In the next section we 
will investigate whether this is a real feature or an 
artifact of our data analysis. Let us concentrate, 
here, on the 3D shape of the Galactic bulge in a 
qualitative way. 

The panel at b = —4° also shows the G alactic 
bar, as traced by iRattenburv et al.l pOOTI ) using 
OGLE II data for stars at a similar latitude. In 
the work of Rattenbury et al. the derived bar was 
arbitrarily shifted in distance so that its center 
would be at 8 kpc. The angle between the struc- 
ture we find and the line of sight is clearly the 
same as that of the Rattebury's bar. The center 
of structure, at 6 = —4°, is '-^ 7 kpc away from the 



Sun (see below). 

In Fig. m each panel shows a vertical section of 
the density map, parallel to ^ = 0° axis, at a given 
longitude. The X — b plane passing from / = 0° 
is in the central, middle panel. The Sun would be 
at {X, b) = (0, 0°), outside each panel, on the left. 
Lines of sight at different latitudes are shown as 
horizontal color strips, in each panel. Again, the 
line of sight at {l,b) — (0°, +4°) is missing due to 
high extinction. Thin white lines in each panel are 
lines of constant height above/below the plane {Z 
linear coordinate). 

The color map has been normalized so that 
the peak density along each line of sight gets the 



maximum yellow intensity. Again, this allows us 
to trace the bulge shape even far away from the 
plane, but of course the density of bulge stars at 
|6| = 4° and |6| — 8° is not the same. Looking at 
\l\ < 1°, the bulge splits in two components both 
above and below the plane. However, when mov- 
ing to sections at positive longitudes (top panels) 
only the one closer to the Sun is visible, while at 
negative longitudes (bottom panels) only the far 
one remains. 

Although we do not have data in the region 
close to the plane, it is reasonable to assume 
that the elongated structures seen at |5| > 4° 
would merge in the center, completing an X- 
shaped structure. If we make the exercise of join- 
ing the arms of the X above and below the plane, 
we find that the whole structure is approximately 
centered at a distance of 7.4 ± 0.6 kpc from the 
Sun, which is consistent with the most recent de- 
termination of the distance to the Galactic center. 
The latter was found to be 7.52 ± 0.10 kpc by 
Nishiyama et al. (2005), 8.0 ± 0.6 kpc by Ghez 
et al. (2008), and 8.33 ± 0.35 kpc by Gillessen et 
al. (2009). The accuracy of the distance from the 
Sun to the center of the X-shaped structure found 
here depends on how appropriate is the adopted 
value for the absolute magnitude of the bulge RC. 
Uncertainties on this parameter may come from 
the broad bulge metallicity distribution function, 
and the unknown precise age and helium content 
of the bulge. However, this is beyond the purpose 
of the present paper, which addresses the shape of 
the outer (|5| > 4°) Galactic bulge. 

In closing this section let us comment on the ef- 
fect of incompleteness. Incompleteness affects the 
present analysis in two ways. First, it artificially 
reduces the star density, especially in the faint RC 
(Fig. [2]). For a qualitative analysis of the bulge 
shape such as the present one, this might not be 
a serious bias. Second, severe incompleteness pro- 
duces a cutoff in the middle of the RC, so that the 
peak in the RC luminosity function looks brighter 
than it really is. Indeed, when fields at |6| < 3.5° 
were tentatively included in the analysis, the map 
would show the bulk of bulge stars getting closer 
and closer to the Sun, for lower and lower |6|. This 
effect was not trusted as real, and this is where our 
decision to exclude fields at |6| < 3.5° comes from. 
Even for high |6|, however, the extinction is highly 
variable on a field to field basis, and so it is the 



2MASS degree of incompleteness. The wiggles of 
the arms of the X in Fig. |4] (see, e.g., the map at 
positive latitudes, for I ~ +5°) might in fact be 
due to differential extinction and thus incomplete- 
ness. 

4. Relative density of the bright and faint 
RC populations 

Here we investigate the density of RC stars in 
the faint and bright RCs. The RC densities not 
only reveal the structure and symmetry of the X- 
shaped bulge, but also provide a consistency check 
to address any lingering doubts that the brightness 
difference between the two RCs might be due to 
something other than distance. 

We measured the density of RC stars as a func- 
tion of vertical height below the Galactic plane, at 
longitude / ~ -|-5° for the bright red clump, and at 
I ^ —3° for the faint red clump; at these latitudes 
the bright and faint RCs, respectively , are max- 
imum fsee iMcWilliam fc Zoccaiill2010l) . The RC 
densities were measured from 30 arc-minute radius 
circular 2 MASS fields, in K^ frequ e ncy h istograms 



similar to iMcWilliam fc Zoccalil (|2010h and 



eluded a linear interpolation to remove the back- 
ground counts. Background subtraction would be 
significantly improved with the use of isochrone 
fitting techniques, but this is beyond the scope 
of the current paper. To compute the space densi- 
ties of the RC stars from the number counts it was 
necessary to assume distances for the bright and 
faint RC populations; we adopted 6.4 and 9.0 kpc 
distances for these populations, respectively (see 
McWiUiam fc Zoccali 2010). 

The vertical extent of our density fields ranges 
from b = —6.5° to —12.6° in latitude, and gener- 
ally corresponds to vertical distances, Z, from 1.0 
to 1.4 kpc. All projection effect have been taking 
into account when computing the stellar surface 
densities for the bright and faint RC populations. 

Our results for the two regions are plotted in 
Fig. [5l which shows that the faint and bright RCs 
possess similar declining trends in space density 
with distance below the Galactic plane. How- 
ever, if one considers only the angular extent below 
the plane (i.e., latitude), it is clear that the faint 
RC density declines much more rapidly with an- 
gle than the bright RC. For example. Fig. [5] shows 
that at 6 = —9° the faint RC has almost vanished. 



but there are large numbers of RC stars in the 
bright RC at this latitude. 

Put in another way: with the assumption of 
approximate foreground-background symmetry of 
the bulge, the trend of RC density with latitude 
is consistent with a distance interpretation for the 
bright and faint red clump populations. Given 
the extensive discussion of distance interpretation 
in IMc William fc Zoccah (2010), and the consis- 
tency check discussed here, we conclude that the 
RC brightness differences must, indeed, be due 
to distance. Therefore, the foreground and back- 
ground bulge RC densities possess approximate 
symmetry. In detail, however, our measured den- 
sities are slightly higher for the faint RC field. We 
believe that this difference is probably real, but 
could be due to errors in our analysis, such as poor 
background subtraction, or in the precise distances 
adopted for the two RCs. Clearly, this small den- 
sity asymmetry should be investigated further. 

5. Discussion and Conclusions 

The position of the RC(s) was mapped across 
a large bulge area, using 2MASS photometry. By 
using the RC(s) magnitude as distance indicator, 
density maps were constructed along the line of 
sight, for the region \l\ < 8.5° and 3.5° < |6| < 
8.5°. 

The main conclusion of the present paper is 
that in the Milky Way bulge has an X-shape, seen 
almost edge on, with an inclination of ^ 20° with 
respect to the line of sight. The X is clearly seen 
within |Z| < 2°, with the closer arm remaining at 
positive latitudes while at negative latitudes only 
the farther arm is present. The arms of the X 
merge together both above and below the plane, 
for latitudes |6| < 4°. Therefore, in Baade's Win- 
dow at 6 = —4° one would just see an elongated 
structure, corresponding to the two lower arms be- 
ing almost completely merged. This feature has 
been classically interpreted as the Galactic bar. 

Our crude measurement of the space densities 
of RC stars as a function of vertical height below 
the Galactic plane is consistent with the adopted 
RC star distances, if the bulge structures studied 
here are roughly symmetric from front to back. 
This removes any lingering doubt about the inter- 
pretation of the RC iiT-band magnitudes as due to 
distance, and strengthens the conclusion that the 
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Fig. 5. — Measured density of RC stars in 
the foreground (black) and background (red) 
RC peaks are shown; units for N^ic are 
10""^ pc~^. Foreground and background RC dis- 
tances of 6.4 kpc and 9.0 kp c were adopted from 
McWilliam fc Zoccalil ()2010l ) , who estimated them 
from b = —8°. Note that the two RC popula- 
tions show very similar declining density trends 
with vertical, Z, distance below the plane, but are 
markedly different when compared with angular 
distance (i.e., latitude) from the Galactic plane. 
In particular, at latitude b = —9° the background 
RC population is present, but barely detectable, 
yet the foreground population at the same lati- 
tude is enormous; only a.t b — —12.6° does the 
foreground RC density resemble the background 
at & = —9°, where Z = 1.425 kpc for both popu- 
lations. This supports our distance interpretation 
for the brightness difference of the two RC popula- 
tions, and indicates that there is approximate den- 
sity symmetry between the foreground and back- 
ground bulge. However, a slightly higher density 
of stars may be present for the background, faint, 
RC. 



bulge contains an X-shape structure. Our densi- 
ties are slightly higher for the background struc- 
ture; this possibility should be investigated with 
more sophisticated techniques. 

Several galaxies are known to host X-shaped 
bulges (e.g., NGC 128, NGC 3625, NGC 4469), 
and some numerical simulations predict their 
formation, with a particul a rly strong buckling 
of a bar (Patsis et al.l 2002: Athanassoulal 200 



Mihos et al. Ill995l ). iMcWilham fc Zoccalil (|2oT 
review some of the proposed mechanisms to form 
such structures. 

The X-shape structure seems to be relatively 
concentrated, with respect to other galaxies (e.g., 
NGC 4710) and probably is the dominant com- 
ponent in the inner 2 kpc. It remains to be un- 
derstood how the bulge 3D structure can be rec- 
onciled with the other features, suc h as the pres- 



ence of a radial metallici ty gradient (jZoccali et al. 



20081 : iBrown et al.ll2010l) , the kinematic difference 
between the metal po or and metal rich compo nent 



in Baade's Window teabusiaux et al.l 120101 ) and 



the behavior of the alpha element ratios, which 
seems to be constant in different fields (Gonzalez 
et al. 2011). It should be kept in mind, however, 
that all the properties listed above have been de- 
rived from observations on ^ ~ 0° only. Just as 
the shape turned out to be much more complex as 
soon as we analyzed a larger area, the other prop- 
erties might turn out to be a partial view of the 
global picture. 

Radial velocities of stars in the two RCs, at 
b = —8° revealed no differences (JDe Propris et al 



20111 ). However a similar analysis at & = —6° re- 
veals the presence of two peaks at Vr = -1-100 / — 
100 kms~^ in the faint /bright RC, respectively 
(Vasquez et al. 2011, in preparation). Proper mo- 
tions also do not seem to differ between stars in 



the two clumps (jVieira et al.ll2007l ). 

Further investigations are clearly necessary to 
completely understand our bulge structure, kine- 
matics, chemistry and, ultimately, origin. The 
present mapping traces the X-shaped structure 
across the whole bulge area and may help to guide 
further analysis, as well as to design future obser- 
vations of the Milky Way bulge. 
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